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PREFACE 


W HILE lecturing in 1913 on the calculation of self and mutual 
induction, the writer noticed that the integral in the case of co- 
axial circles was of such a form as to make the application of the 
arithmetico-geometrical means obvious. An extremely elegant formula 
{Proc, Boy, Soc. A, Vol. c. 1921, p. 63. Appendix, this book, No. 35) 
when put to numerical test justified the method as one well adapted to 
computation. The calculation of more complex cases of self and mutual 
induction led to several interesting formulae for the complete third 
elliptic integral. With a view to making a complete study of the group 
of recurrence formulae associated with the a.g.m. scales, further work 
was postponed until a vacation in 1921, when many of the results in- 
cluded in this book were obtained, and were thought by the writer to 
be new. On later consulting the Collected Works of Legendre, Gauss 
and Jacobi, the writer found that he had rediscovered many known 
formulae, but that quite a number, likely to be of service in computation, 
appeared to be hitherto unknown. In the summer of 1922 the writer 
decided, as a vacation task, to include in a single monograph the entire 
theory of elliptic functions associated with the a.g.m. scales, thus com- 
pleting, and to some extent adding to, the work of Gauss left unpublished 
after his death, and placing before mathematicians in accessible form 
a mode of approach to Elliptic Function Theory directly related to the 
art of machine-computation. It is interesting to note in this connection 
that this subject was to have formed the content of the third volume 
of Halphen’s Fo^ictions ElliptiqueSy unfortunately left incomplete at the 
time of his death. 

While the present monograph contains many new formulae and 
methods of computation, no claim is made as to novelty in fundamental 
analytical treatment. References to more important books and memoirs 
consulted by the author are given in foot-notes. A complete bibliography 
is not given since the reader will find exhaustive references in the Boyal 
Society Index, 1800-1900, Vol. I, Mathematics, and after 1900 in the 
mathematical volumes of the International Catalogue of Scientific 
Literature, 
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ON THE DIRECT NUMERICAL CALCULATION 
OF ELLIPTIC FUNCTIONS AND INTEGRALS 

SECTION I 

INTRODUCTION 

Although the numerical evaluation of elliptic functions and integrals 
is dealt with in many of the standard treatises, it cannot be said that 
the use of tables is entirely satisfactory in dealing with computations 
required in several branches of physics and astronomy. Many of the 
integrals appearing in these problems can only be expressed as some- 
what complicated expressions involving the complete elliptic integrals 
of the first and second kinds, and the corresponding functions of com- 
plementary modulus. These formulae give rise in many cases to 
differences of nearly equal quantities : the numerical calculation of such 
expressions to a sufficient number of significant figures is thus often a 
matter of some difficulty, requiring tedious and sometimes uncertain 
interpolations. 

While the numerical calculation of the Theta-functions is extremely 
rapid by the use of the highly convergent series involving the name 
these are not usually the functions which make their appearance directly 
in physical or astronomical problems. The necessary transformations 
required to express the results in terms of the g-series require in many 
cases somewhat complicated analysis, while numerical computation by 
this method necessitates the exact calculation of a large number of 
auxiliary quantities*. 

The forthcoming tables in course of construction under the auspices 
of the British Association t are so arranged as to allow a large number 
of functions to be evaluated, and will very greatly aid numerical work 
in many practical problems involving the use of elliptic integrals and 
functions. While sufficient for most purposes, the use of these tables 
will still involve considerable interpolation when it is required to take 
out values corresponding to arguments and moduli expressed by six or 

* On this point see Kosa, E. H., and Grover, F. N., Bulletin 169, Bureau of Stan- 
dards, Washington (1913), pp. 67 and 73; Nagaoka, H., Tok. ColL Sc, J. 27 (1909), 
No. 6, and Jour, Tok. Math, Fhys. Soc. 4 (1908), p. 284; 6 (1911), p. 10; Nagaoka and 
Sakurai, Scientific Papers of the Institute of Physical and Chemical Research, Vol. ii. 
Dec, 1922, pp. 1-67; Olshausen, G. R., Phys. Rev. 51 (Dec. 1910), pp. 617-636. 

t See report (by Sir A. G. Greenhill), Brit, Ass. Report, 1912 (Dundee), pp. 
39-56. 
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seven significant figures, as is now necessary in several modern applica- 
tions of elliptic functions. 

For several reasons it is advantageous to have available a direct 
method of computation, independent of auxiliary tables, based on the 
use of the modern high-speed calculating machine. The writer has 
therefore thought it desirable to set out in detail the method of 
calculation based on Landen’s quadric transformation and the use of 
the scale of arithmetico-geometrical means, recapitulating briefly the 
original developments of Lagrange, Legendre, Grauss, Jacobi and others, 
extending them in several directions and keeping in mind the ultimate 
use of the calculating machine for numerical work. 


SECTION II 

HISTORICAL NOTE ON LANDEN^S TRANSFORMATION AND 
THE VARIOUS SCALES OF MODULI AND AMPLITUDES 

The importance of realizing rapid and accurate methods of calculating 
the elliptic integrals now denoted by 

u = ^)= f (^, k) and E{4>y /fc) = T A (<#>, k) . . .(1), 

jo 

where A (</>, h) sin^ <#»), 

was first remarked by Euler* (1766), although it was not until several 
years later that Landenf (1775) discovered in geometrical form the 
transformation which forms the basis of the existing methods of the 
numerical calculation of the elliptic integrals. 

A method of successive transformations for the ultimate reduction of 
the algebraic forms of these integrals to elementary integrals was pub- 
lished by Lagrange t in 1784-85. This memoir contains an exposition of 
the scales of arithmetico-geometrical means, together with two types of 
algebraic transformations corresponding to the increasing and diminishing 
amplitudes of Landen’s trigonometrical forms: it is shown that the 
given integral must lie between two limits which may be made to 
approach each other as closely as may be desired : the limits thus ob- 

* Euler, Novi Comm. Acad. Sc. Petrop. vol. x. 1766, 

t Landen, Phil. Trans. Hoy. Soc. vol. lxv. p. 283, 1775; Mathematical Memoirs, 
London, 1780. 

+ Lagrange, “Sur une nouvelle mdthode de calcul integral pour les diff^rentielles 
affect^es d’un radical carrd sous lequel la variable ne passe pas le quatri^me degrd,” 
M6m. de VAcad.roy. des Sc. de Turin, t. n. 1784-5; (Euvres (Gauthier- Villars, Paris, 
1868), t. II. pp. 253-312. 
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tained for the integrals of the first and second kinds are practically 
identical with the limits and series obtained about the same time by 
Legendre. 

Landen's transformation was applied by Legendre* to the numerical 
calculation of the elliptic integrals, and furnished the method of com- 
putation by means of which the latter's well-known tables were con- 
structed. The scales of moduli and amplitudes employed by Legendre 
are briefly described below. 

If we form successive moduli and amplitudes according to the 
recurrence formulae 

K^i=^{l-kn)l{l+Jcn) and tan(0,,+i-<^,,) = ^,,' tanc^,, ...(2), 
Landen’s transformation in trigonometrical form leads to the result 

= (1 + JCn) K) ( 3 ). 

As n increases, the successive moduli kn converge to zero, and the 
complementary moduli kn to unity, while the amplitudes increase so 
that tends rapidly to a definite limit : ultimately we may write 
F(cl)n, kn) - <i>n and thus obtain the result 

h) = (2A'/7r) <l>nl2\ where Z’= Jtt (1 + ^i) (1 4 - • • • (1 + ^n) 

(4). 

The result of calculating successive moduli and amplitudes in reverse 
order is equivalent to calculating in forward order these quantities 
according to the recurrence formulae 

kn+i=2jknl(l-^kn) and sm(2ll/n+i-\l/n)=kn&m^n •••(5). 

This transformation then gives 

F(}l/n+ii kn+i) — J (1 + kn) F(\l/n, kn) (6). 

As n increases the successive moduli kn+i converge to unity and the 
amplitudes to a limit if/n' we thus obtain, when n is suflSciently large, 

F(}l/n,kn)= I seci/^n(^ifrn = logtan(i7r + ^if^n) ( 7 ), 

giving finally 

F k^ = sj (kik^^z • • • kri)! sf k ^ . log tan -f ^^n) • • • (3)* 
Legendre also developed for purposes of numerical calculation extremely 
convergent series for E(i>, k) which he employed in the tabulation of 
this function. 

* Legendre, “M^moire sur les integrations par arcs d’ellipse” and “Second 
memoire, etc.,” M^m, de VAcad. des Sciences de PariSy ann. 1786 (Paris, 1788), 
pp. 616-643 and 644-683 ; TraiU des Fonctions Elliptiquesy Paris, 1825, t. i. 
p. 79 et seq. 
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The method discovered independently by Gauss and called by him 
the “algorithm of the arithmetico-geometrical mean'' originated in 
connection with the evaluation of complete elliptic integrals arising 
from a problem in attractions required in planetary theory*: in a note he 
mentions that the results were obtained by him as part of a more com- 
prehensive theory, independently of the results of Lagrange and Legendre 
with which they are closely connected. Gauss employed a trigono- 
metrical transformation along the following lines. 

Commencing with two numbers {a^, successive numbers hri) 
are calculated from the recurrence formulae, 

^ ^m.+1 “ fj (^n^n) (^)* 

Writing cos^ + 6/ sin^ Gauss employs the recurrence 

formula 

an tan tan <t>ny 

or tan = tan + K)/ ^n) •••(lO;. 

From the relation 

tan (t>n = (AjAg ... A^. ..)/(aoai . .. . ..) . tan 

we have ultimately, since the <^'s tend to a limit <^n, 

r<#>o 

Jo 

and hence 

a„ = tan-* [( Aj A, . . . A„ . . . . . . ) • tan ... (1 1). 

Extremely convergent series for the elliptic integral of the second 
kind were also obtained by Gauss. The application of this method to 
the numerical calculation of the elliptic integrals was considered by 
Jacobi some years later +. 

In the course of researches carried out between 1797 and 1808, but 
not published until after his death considerably later Gauss developed 
the theory of elliptic functions to a remarkable extent and derived a 
large number of formulae which, at a later date, were obtained inde- 
pendently by Abel and Jacobi and which are now associated with the 

* Gauss, “ Determinatio Attractionis quam in Punctum quodvis Positionis datae 
exerceret Planeta, etc.,” Comm, Gott, Soc. Beg, Scient, iv. 1818; Oes, Werke 
(Gottingen, 1866), Bd. in. pp. 331-355. The note mentioned above is dated Feb. 9, 
1818 ; Ges. Werke, Bd. iii. pp. 357-360. 

t Jacobi, Fundamenta Nova, 1829, sections 38 and 52 ; Ges, Werke (Berlin, 1881), 
Bd. I. pp. 154 and 203. Also “Nnmerische Berechnung der Elliptisohen Functionen,” 
Crelle, Bd. xxvi. pp. 93-114; Ges. Werke (1861), Bd. i. pp. 345-368. 

X Gauss, Ges. Werke (1866), Bd. in. pp. 361-403 and 433-469. 
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theory of the Theta-functions. In one of these posthumous papers, Gauss* 
employs the recurrence formula tan2</>n= \/(^n+i/a«+i) tan which 
is much more convenient for purposes of numerical calculation than (10). 
The last mentioned method is considerably simpler than that of 
Legendre’s scales of moduli, especially when a calculating machine is 
available, although in theory these two methods of computation are 
practically identical. 

Although a few formulae based on the scales of arithmetico-geometrical 
means were briefly considered by Jacobi t, no attempt seems to have been 
made by him to apply the results to actual computation. 

As the method of evaluating the elliptic functions by the use of the 
A.G.M. scales is only briefly dealt with in existing text-books J, and does 
not seem to have been developed to the extent it deserves in various 
memoirs which have appeared on the subject in recent years, a tolerably 
full account of the subject is given in the present book, together with 
a brief outline of the main analytical steps. Several new formulae have 
been obtained by the writer which facilitate the application of this 
method of numerical computation and considerably extend its scope. 

From the point of view of pure analysis, the modern Weierstrassian 
notation presents many advantages in regard to elegance and symmetry, 
a superiority which no longer holds, however, when it is required to 
proceed to numerical evaluation. As the main point of the present 
paper is concerned with improvements in methods of computation, the 
older Jacobian notation has been adhered to as far as possible, especially 
as it expresses most readily the elliptic functions and integrals which 
arise from integrations involving circular and hyperbolic functions. The 
results are set out in such a manner as to provide a fairly complete 
compendium of formulae intended to be of service to the computer in 
reducing to numerical results such expressions involving elliptic functions 
and integrals as may arise in physical and astronomical problems. 

* Gauss, loc. cit. p. 388. 

t Jacobi, Oes. Werke (1881), p. 357. 

Z Cayley, A., Elliptic Functions, 2nd Ed. (G. Bell and Sons, London, 1896), 
pp. 326-338; Enneper-Miiller, Elliptische Functionen (2nd ed. 1890), pp. 361-364; 
Tannery-Molk, Fonctiom Elliptiques, t. iv. 1902, pp. 269-276. According to a foot- 
note in Halphen*s Fonctions Elliptiques, t. ii. 1888, p. 310, the subject of the 
arithmetico-geometrical mean was to have received special attention in the author’s 
third volume, unfortunately incomplete at the time of his death. 
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SECTION III 


ON THE SCALE OF ARITHMETICO-GEOMETKICAL MEANS 


In forming the scale of arithmetico-geometrical means, we start with 
the positive numbers (uo» b^y of which ao is the greater, and form 

successively 
ai = K«o + &o), 

6i = V(«o&o), c, = |(ao-6o) 1 


Oa = K«i + 

bi = J{(hh), = 

1 

— J(%-1 "b bn^i) 

bn ~ ^Jif^n-^lbn-T^ \ {P'n-i “ 5ot,_i) j 

i 

We obviously have the relations 



«» = «^«+i + c„+i , = a,^ - b/\ 

bn~ ^n+l~ ^n+li 

(13). 

The a's and t’s tend to the same limit denoted by 



M(aoybo) - lim «« 

(14) 




with extraordinary rapidity, even when ao and bo are initially numbers 
of very different magnitudes*. 

We notice in passing that if (^o, bo) satisfies a homogeneity relation 
of the form 

€M(aoi bo)=^M(faoy ebo) (15), 

€ being any number. The array of numbers (12) will be referred to as 
the ‘‘ scale of arithmetico-geometrical means (^o, or more briefly 
as the A.o.M. scale bo). Also, the limit if (ao» h) will be denoted 
by an, as long as by so doing no ambiguity is involved. 

If we calculate in the same way an array of numbers commencing 
with {ao =aoy bo -Cq- fJ(ao -'bo)}y we derive what will be called the 
'‘complementary a.g.m. scale {a^, 5o'),” the symbols being in this case 
denoted by accented letters, and the limit Micio, bo) by a^. 

* The usual arithmetical process of extracting square roots is easily adapted to 
most calculating machines of modern type. An account of these is given in 
Modem Instruments and Methods of CalculatioHy edited by E. M. Horseburgh 
{London, G. Bell and Sons, 1914) ; also by d’Ocagne, Le Calcul Simplifii (Gauthier- 
Villars, Paris, 1905). In carrying out the process the work is considerably shortened 
by keeping in mind the rule that if the first p digits out of the number n 
required in the square root have been obtained by the usual process, the next 
p-1 digits can be obtained by division only, with a possible error of 1 in the last digit 
(ChrystaPs AlgehrOy Part i. 5th ed. 1904, p. 210). To illustrate the rapid convergence 
of the A.O.M. scales. Gauss works out an example starting with ao=l and 6o=0*2 : 
ag and bg are represented by the same number to fifteen significant figures. 
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It is obvious that the operations of the a.g.m. scale may be carried 
out backwards : if the entries are denoted by negative suffixes, they 
are connected with those of the complementary scale by the relations 

— — (l^)’ 

These relations evidently hold good when accented and unaccented 
symbols are interchanged. As a result of these relations it will be 
seen that there is no advantage to be gained from considering any other 
than the a.g.m. scale of positive suffixes (an, bn) and the corresponding 
complementary scale (a^ , bn)- 


SECTION IV 

LANDEN^S SCALE OF INCREASING AMPLITUDES 

(i) Calculation of the Elliptic Integrals of the First and Second 
kinds, modulus k. 

If we write A,j= v^(an^cos^<^n + ^«^sin^<j!>n), the recurrence formula 


tan (<^n+l - <#»n) = (bnldn) tail <^n *(17) 

may be written in either of the forms 
sin (2<^n - sin 

COS(2<^n~<^n+l) = An+l/«n+l (lb). 

From these we derive 


An+I + Cn+i COS <^n+i = and An+1 - Cn+1 COS <#)n+i = dnbjK - • -(19), 
which give on differentiation, 

#o/Ao = i#i/Ai = |d:<^2/A2= ... = (J)^#n/An= (20). 

From the relation 

AVl + Cn+iAn+i COS <^n+l = ^ An^ + ^Gubn 

we obtain, making use of (20), 

An<^^n “ ^n-¥ld<hn+l ~ ^n+i COS ^n+ld<hn+i ^nl^nd4*n / • • •(21)* 
Subtracting from each side the identity 

a>f^d<hnl^n Of^n+ld^n+ij ^n+1 ~ d<f>nf^n ■" dnbnd^pnl • • *(22) 
we obtain 

(An anl^i^d<hn (An+i ”* a^n+ij ^n+i^d<f>n+l 

— Cn+i COS ^n+l^^n+i ” Wd<l>nlK ...(23). 
Writing down a series of such equations, commencing with n = 0, 
and noting that as n increases, (An- Unj^n) tends to zero in the limit, 
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we have, on integrating, 

f<l>0 

I Aorf<5f>o~ao® / c?<^o/^o = Cisin<^i + C2sin<^2+ ... +c,isin </>„+ ... 

Jo Jo 

- i (co” + 2c, + 4c/ + . . . + 2’*c/ + . . (24). 

If we now construct the a.g.m. scale (ao=l, bQ = k'), and calculate 
successive amplitudes commencing with = we have Ao = A(0, ^). 
As n increases, the a's and b's tend to the same limit while the 
angle <bnl2^ tends to a finite limit. We thus obtain from (20) and (24), 


u^F{<i>,k)=^(llan){<i>nir) (25), 

Jc)-F(Ji, k) = {ci sin + C2sin^2+ ... + sin + ...) 

+ Aci + ... + 2^c^ + ...)F(<l>y k) (26). 

If we commence with <^o = we have 

</>2 27 r , <^3 = 47 r ... = Jtt , 

SO that K = and 


(K-’ E)/K = + 4c. 2^ + . . . + 2^c,,2 + , , (27). 

The series (24) and (27) converge with extreme rapidity, resembling 
in this respect the g-series with which they are closely connected*. 


Jacobi’s integral Z (<b, k) is defined by the relation 

Z(<b, k) = k) ^ (EIJT) F{<t>, k) (28) 

or, if the angle <#» is connected with the argument u by the relation 

sin <j(> = sn (^^, k), 

Z{u,k) = (u, k) du - (E/K) u (29). 

From (26) and (27) we derive the expansion 

^(<^, ^)=Cisin<^i + C2sin<?^2+ ... +c„sin<^„ -f (30). 

In terms of Jacobi’s Theta-function ©(^^, k), 

Z (u, k)du = Z(<^, k) d<l>l^ (<^, k) = log [@ (u, k)/(^ (0, k)] . . .(31) 

where 0 (0, k) = J(2k'Kl7r) (32). 


* These relations are discussed by Legendre {TraitS des Fonctions ElliptiqueSy 
t. in. Deuxi^me Supplement, 1828, p. Ill et seq.). On the connection between the 
A.o.M. scale and the modern transformation theory of the Theta-functions, see 
Tannery and Molk, Fonctions Elliptiques, t. iv. 1902, note 2, pp. 269-273. A develop- 
ment of the theory of the Theta-functions based on Legendre’s scales of moduli is 
outlined by Richelot, Correspondenz mit Herrn Professor Sehrbter (Ednigsberg, 
1868). See Examples 19-24 of the Appendix. 
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Prom (20) we have = so that the series (30) 

enables us to write 

f Z k) di>l^ (<^, ^) = 2 (Pn/2^) f ” sin 
•'O 1 Jo 

On integrating and making use of (19) and (18), we obtain 

^ sin €f>nd<f}n/K = (l/Cn) log {an-ij ^n-^i) = {ijCn) log SCO (2</>n-2 ~ S^>n-l), 

leading to the expansions* 

( Z{Uy k) du ~ 2 (l")^ leg («»-i/^n-l) (^^)> 

/o 1 

or 

log [® (u, k)l% (0, ky\ = J log (^o/^o) + i log sec (20^ - ^i) 

+ 1 log sec (2<^i - </)2) + . . . + (^)”+i log sec (2<^,,.i - <^n) + (34). 

In the series (34) it is more convenient to write 
ajA^ = cos (<^x - <^o)/cos 00 
for purposes of numerical calculation. 

(ii) Calculation of the Jacobian Functions sn k), cn (w, A), dn (w, A:) 
etc. in terms of the argument u. 

When the argument u is given and the modulus k is known, it is 
only necessary to compute the a.g.m. scale = ^o such a 

value of n that Cn is less than the small quantity determining the order 
of accuracy of the calculations. Since u = F{(f» k), equation (25) enables 
us to calculate in circular measure the angle (hn from the relation 

<hn^T^anU (35). 

From the recurrence formula (18), 

sin (20n_i - 0n) = (Cnian) siu 0n, 

we are enabled to calculate successively the angles <hn-i •••02, 0i and 
finally 0o, in terms of which 

sn (uy k) = sin , cn {% k) = cos 0o, dn {u, k) = Aq - cos 0o/ cos (0i - 0o) 

(36). 

Z {% k) may then be computed from the series (30), and @ {Uy k) from 
formulae (32) and (34). 

^ Formulae (30), (33) and (34) are practically the same as those given by Legendre 
{Traits des Fonctions Elliptiquesy t. i. section 90), and are contained implicitly in 
formulae given in Gauss’ posthumous papers. Formula (34) is equivalent to that 
given by Jacobi (CrelUy Bd. xxvi. pp. 93-114; Gesammelte Werke (1881), Bd. ii. 
p. 357). 
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SECTION V 

THE HYPERBOLIC SCALE OF INCREASING AMPLITUDES 

(i) Calculation of the Elliptic Integrals of the First and Second 
hinds^ modulus 

Writing i^hn for throughout the formulae (17) to (24), we obtain, 
instead of (17), 

tanh (<^n+i - <l>n) = (K/Un) tanh <l>n (37), 

which may also be written 

sinh (2<l>n - ^n-fi) = (c»+i/«n+i) sinh (hn+i , cosh (2<^„ - = ^n+l/(^n+i 

(38). 

If we denote lin= Ay(a„^cosh^<^«--i„“sinh^<^,i), it is easily proved, as 
before, that 

^<^*^ 2 /^ 2 = ••• —(^}frd<j>nl^ii^ (39) 

from which it follows that as n increases, we have in the limit, 

[^'^c^o/Ao = (1K)(<^J2-). 

Jo 

If we now write 

sinh = tan (40) 

it is easily proved that 

Ao = A (<^, ^')/co8 <l>, rf<^o/Ao = d^/A (41), 

SO that 

u^F{<hy h:') = (IK) (<^«/2”) (42). 

Formula (24) becomes, 
r0o - 

I ^Qdi >^ - / d<t>jA^ - (cj^ sinh <l>i + sinh <^2 + • • • + sinh <l>n+ 

Jo Jo 

r^o 

-|(c„“ + 2c/+4c/+.-.+2’*c„“+...) ...(43). 

Jo 

It is easily shown that 

f ' { A(4>,t)d4>fcos!‘<f> 

Jo Jo 

ru ru 

= / dn® (w, k') dufon^ (u, k‘) = / dn® (m, k) du, 
Jo Jo 

or in terms of the ^-function of imaginary argument, 

I Aq d<ho {1 1 i)Z {iu, k) + (EjK) u, 

Jo 

Substituting in (43), and making use of (27), we derive the series, 
corresponding to that of (30), 

(l/^)i2^(m, ^) = Cl sinh ^1 + 02 sinh 02+ ••• + sinh ... ...(44). 
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The following formula, corresponding to (34), easily follows by 
integration, 

log [© (^^4, A)/@ (0, k)] = i log (ajA^) + 2 log sech - i>n) 

(45) 

in which we may conveniently write aJA^ = cosh (<#>i - <^o)/cosh <#>o. 
Prom (42) we easily deduce the relations, 

(l/i) sn (m, k) = sn (u, A')/cn (u, K) = sinh <^o 1 

cn (m, >J:) = 1 /cn (Uy k') = cosh <^o }■ 

dn (iuy k) ~ dn (u, k')lm (Uy k') = Aq = cosh <^o/cosh (<^i - <j!>o) J 

(46). 

Making use of Jacobi's well-known relation in the form 

Z (Uy k') = sn (uy k*) dn (^^, i')/cn (w, A') - iruK^KK') + (^^^, A) 

(47) 

we have from (28) and (46) 

AO 

= sinh <?!>o/cosh (<^i - <^o) - '^F(<l>y K)I{2KK) + iZ{iUy k). 
Noting that Legendre's relation gives 

EIK + EIK'~l = \itl{KK) (48) 

we obtain from (27) 

E{<f}y k') = sinh <^o/cosh (<^i — <^o) + (1 — EjK) F{<t>y k') + iZ (iu, k)y 
or, finally, making use of (27), (42) and (44), 

E{<l>y k') = sinh ^o/cosh {<l>i - 

+ |(co^ -h 2c/+ ... + ...) (IK) 

— (cj sinh cpi + C 2 sinh ^2 + • • • + sinh <^„ + . . .) ... (49). 

(ii) Calculation of the Jacobian Functions sn {Uy k'), cn (w, k')y 
dn (Uy k')y etc, in terTns of the argument u. 

Compute the a . g . m . scale {a^-ly bo = k') and calculate <pn from the 
relation being such that Cn is less than the small quantity 

determining the order of accuracy of the calculations. The hyperbolic 
recurrence formula 

sinh = (cnK) sinh 

enables us to calculate successively ... ^ 2 , and finally <^o> in 
terms of which we have from (46) 

m{Uyk') = tmh.(pQy cn{Uy k') = l/cosh<l>Qy dn(t«, i0 = 1/cosh (^i-~<^o) 

(50). 

Z{Uy K) may then be calculated from (47) and (44). 
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(iii) Calculation of the Jacobian Functions of imaginary argument^ 
sn h\ cn {iu, A:), dn (m, K), etc. 

The procedure outlined above, together with equations (46), obviously 
enables us to calculate the Jacobian functions of imaginary argument 
and modulus 4 while formulae (44) and (45) serve for the functions 
Z (iuy k) and @ (m, k). 


SECTION VI 

GAUSS’ SCALE OF INCREASING AMPLITUDES 

(i) Circular Scale. 

In connection with the a.g.m. scale, Gauss made use of the recurrence 
formula* 

tan Xn+i = sJiO'n+ilK+i) tan 2x« (51). 

It is easily proved that the above formula is connected with Landen’s 
scale (17) by the relation 

tan Xn = sjibnjan) tan (52) 

or 

^n+i — ~ (53). 

In the limit, as n tends to infinity, (bjan) approaches unity, so that 
Xn and both tend to the same limit. It thus follows from (25) that 
if we form the a.g.m. scale (ao = l, bo = k'), and starting with xo deter- 
mined by tanxo= n/(V<^o) tan<^, compute successively Xu X2> ••• Xn ••• > 
we have 

i^(<#>,^) = (lK)(Xn/2”) (54). 

From (52) and (53) it is a simple matter to derive successively 
c^i, ^2 • • • • • • for use in the formulae of Section IV. When these angles 

are required, it is, however, simpler to make direct use of Landen’s 
recurrence formula (17) or (18). 

(ii) Hyperbolic Scale. 

It is evident that formulae (51) and (52) may be written in the 


h 3 rperbolic form 

tanhxn+i = N/(«n+i/^n+i) tanh 2xn (55), 

tanh Xn = fipnlan) tanh <l>n (56). 


* GauBS, Ge». Werke (1866), Bd. in. p. 388. To compare Gauss’ results with 
those of the present paper we note that It then appears 

(p. 392) that formula (26) of the present book was also obtained by Gauss, inde- 
pendently of Legendre (foot-note on p. 9). 
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As n increases indefinitely, and tend to the same limit, so that 
if we write 

tan = sinh <^o or sin = tanh (57) 

it follows from (42) that 

i^(0, F) = (IK) (c#>n/2") = (IK) ixnjT) (58). 

Noting that, as before, 

^n+l~ ^n~^Xn-i .....(59) 

it is a simple matter to derive successively <^ 2 , ••• ••• for use in 


the formulae of Section V. In such cases, however, it is simpler to 
make direct use of the hyperbolic recurrence formulae (37) or (38). 

(iii) Calculation of K\ E' and the name 

If we write <^ = l7r, it follows from (57) that tanh <^o=l, and from 
(56) that tanh Xo = \/(^oK)* From (55) we deduce that 

tanh 2 xo = 6i/ai, tanh2xi = ^ 2 K, •“ tanh 2xri = ^n+iK+i, 

and thus 2xn = i log {(«n+l + *n+l)/K+l - K^i)} = I log {an^^hn^^) 
or (Xn/2’^) = log (an+ 2 lCn+ 2 )- 

It follows from (58) that as n is indefinitely increased, ( J)^ log (ajcn) 
tends to a finite limit given by 

k') = J log { 1 /( 2 ) = {^)^ log {ajcn) . . .(60). 

Making use of the relation (13), anCn = lc\_iy we have 
log (an/Cn) = 2 log {2an/an.i) + 2 log {an~i/Cn,-j)y 
which gives by successive application, 
log(«,iK)- 

2 log {2aJan-0 + 4 log {^a^^^lan-^) + • . . + 2” log {2a,la^) + 2^ log {a^/c^). 

Inverting the order of the series, we obtain 

(iT log (an/Cn) = log (Uo/Co) + log2{l + i + J+ ... + (i)^-^} 

+ log (^iK) + i log {aM) + • . • + (^)’^~' log (ajan^i). 
As n is increased indefinitely we obtain the series 

a.nK' = 1 log (4ai/ci) - 2 (I)’* log (c^^K+i) .(61). 


Making use of Legendre’s relation (48) in the form 

E^.{l-EIK)K' + an 

and using (27), we obtain the following useful formulae, 

E' = \ (l/«n) {c^^ + 2 c/ 4 - 4c2^ + . . . + 2^6-/ + . . .) 

X {log (4ai/ci) - log (o^iK) ~ i log (c^-iK) -■..} + (62), 

r (IK) (1 -W- Cl - 2 c/ - ... - 2 "-^c/ - ...) 

X {log (4ai/ci) - log (^iK) - 1 log (« 2 /%) (63) 
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for the calculation of the complete elliptic integrals of complementary 
modulus in terms of the a.g.m. scale (ao = 1, = ^0- 

The nome q is defined by the relation log (1/g) = 7rjr7if, from which 
it follows that 

log ^ ^ ”■ log (4ai/ci) (^4) 

which is more convenient for numerical computation than (60). It may 
be noted that the entire series of calculations may be carried out with 
ordinary logarithms to the base 10, enabling q and powers of q to be 
evaluated directly*. 

The complementary nome q' = jg obviously related to q by the 

formula 

log q . log q' = TT^. 

SECTION VII 

LANDEN'S SCALE OF DECKEASING AMPLITUDES 

(i) Calculation of the Elliptic Integrals of the First and Second kinds y 
modulus k\ 

If we denote 

+ (65), 

the recurrence formula 

tan (il/n - if'n+l) = (Cn-i-ilan+i) tan if/n+i (66) 

may be written in either of the forms 

= cos = Jan ...(67). 

From these we deduce 

l^n cos if/n — ^^n+li l^n COS ij/n^ 2an+iCn+il^n+l ...( 68 ), 
which give by difierentiation 

d\l/o/\ = #i/Ai = dfj^^ = . . . = dfn/^n = (69). 

If we construct the a.g.m. scale (ao=l, h-k') we notice that as n 
increases, fn tends rapidly to a finite limit. If we start with = 

= A (if/y k')y we have, since c^ ultimately tends to zero, 

* The limit derived in (60) does not appear to be generally known, although an 
equivalent formula is given by Legendre (Fonctions Elliptiques, t. i. pp. 94, 101). 
Although formula (60) does not appear to have been explicitly stated by Gauss, 
the series (61) is found among a number of formulae given in his posthumous 
papers (Gauss, Ges. Werhcy Bd. in. p. 377). Neither (60) nor (61) appears to have 
been noticed by Jacobi or by subsequent writers on elliptic functions. 
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M = ^0 = (l/»n) #n/cos «/»„ = (!/«„) log tan + ...(70), 

a useful formula due to Legendre*. 

From (68) we easily obtain the relation 
+ bn^n COS 

which gives 

~ ^»COS ^n' 

Writing down the system of equations for w = 0, 1, 2, ... w, multiplying 
successively by 1, 2, 4, ... 2^, adding and integrating, we obtain the 
result 


f'I'O , fYn+l r\i/o 

Jo ^odif/Q — 2’*''’'^ ^7i+l<^^n+i = I'CCo^ + 2 Ci^H- 4C2^ + . . . + 2”Cn^) J d\l/o/ A Q 

- (bo sin ilfo + 2bi sin + . .. + sin (71). 

Adding to each side the identity 

/•’/'n+l 

/ ^w+i OOS ^n-hld^n+1 ~ ^n+l siu if^n+i + (2”^"*^^ “^ ^) ^n+i siu t/'n+i 

and re-arranging the terms, (71) may be written 
/•’Ao r^n+l 

Jo Jq (^n+l'~ bn^-iCOSl/^ni-i) dif^n-t-1 

= ^(co^+ 2Ci^+ ... 4-2”c7t^) Jo #o/^o + ^n+isiniAn+i + (^>isini/^i-6oSiniAo) 

+ 3 (bs sin i /^2 - bi sin i/^j) + . . . + (2’''*'^ - 1) (bn^i sin - bn sin i^n) 

(72). 

It is not difficult to establish the identity 
bn+i sin lAw+l -- bn sin f „ = 2Cn+<2 tan (2l/^n-h2 - »/^n+l)/C0S (2i/^n+3 - ^n+ 2 ) 

;(73). 

As 9z increases indefinitely, a^-i-i and bn+i tend to the same limit an, 
while Cn+i tends to zero. It follows that the second integral on the left- 
hand side of (72) tends to zero, so that we finally obtain, making use 
of (70), 

k') = an sin il/n + (Co^ + 2ci^ + ... + 2"c«^ + . ..) 

X (Ijan) log tan (^^ + ^fn) + S (l/^n, On) (74), 

where for the sake of brevity we denote 

^Krnt n) ^ COS (2i/' 3 — i/za) * COS (2i/f4 — l/cs) 


rfn+l 


+ 2(2»-1)c„^3 


tan ( 2 i/'„+ 2 — i/'ra+i) 


‘ COS (2fn^.a- 4'n+O'^ 

Legendre, Fonctiom Elliptiques , 1. 1. chap. xix. 
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a useful formula for the numerical calculation of the second elliptic 
integral*. 

Making use of (27), (28), (48) and (70), we obtain from (74) the 

following interesting expressions 

Z {\f/, K) + = an sin -f S (jA„, Cn) . . .(76), 

2 

or, in terms of u, by (70) we have 


and 


sin xj/n = tanh (anu), | 

Z (u, K) + anujK' = an tanh {anu) + S {\pn, c,,) J 


...(77). 


Jacobi's relation (47), taken with equation (76), gives 

(lli)Z(iu, k) = tan cos (2i/rj - xj/^) - an sin - 2 (iA«, Cn) (78). 

2 

The addition theorem 

Z(u) + Z(v) -Z(u + v)== k^sii (u + v)snumv (79), 

with iu for u and = ± {K + iK') gives two interesting formulae 

iZ(iu + K+ iK\ k) ==an{l + sin lA^) + 2 (lAn, Cn) 

2 

(1/i) Z(iu-K- iK\ k) = an (1 - sin xjjn) Cn) 

2 

If we notice that 



/''Pn 

I sin </'„#„/ = (l/b„) log (a„+i/A„^i) = (1/bn) log sec ( 2</'„+2 - 
Jo 

the result of integrating (72) written in the form 

oo 

Z(u,k') + anU/K' = an sini/'„ + 2 (bn+i sin - bn sin fn) 

0 

gives rise to a formula for the numerical calculation of ® (u, k'\ 

-logsec*.. I (81), 

in which, if desirable, we may write sec xpn = cosh (an^^). 


(ii) Calculation of the Jacobian Functions sn {u, ^'), cn {u^ dn (uy ^'), 

etc. in terms of the argument u. 

Form the a.g.m. scale (ao = l, hf^ — kJ') and determine ipn from 
sin xpn = tanh (anu), n being such that Cn is less than the small quantity 
determining the order of accuracy of the calculations. 

* The series (71) was obtained by Legendre (Fonctions ElUptiques, 1. 1 . chap, xxi.) 
and left by him in this non-convergent form. As far as the writer is aware, the 
convergent series (75) is new. 



LANDKN'S SCALE OF DECBEASING AMPLITUDES 


17 


From the recurrence formula (66), tan(^„_i-^„) = (c„/a»)tan^„, 
compute ... 1 ^ 1 , and finally in terms of which 

sn («, Jc') = sin i/',,, cn (u, k') = cos i/'„ , dn (m, A') = A, = cos (2i/'i - ^o) 

( 82 ); 

(m, k') may then be calculated from (77) and ® (m, k‘) from (81). 


SECTION VIII 

THE HYPERBOLIC SCALE OF DECREASING AMPLITUDES 

(i) Calculation of the Elliptic Integrals of the First and Second kinds, 
modulus k. 

If we write iti/ for tj/ throughout the formulae (65) to (69), (65) becomes 

A» = V(«n‘* cosh’’ fn - Cn^ sinh“ i/f„) (83), 

while the recurrence formula (66) may he written 

tanh (ij/n - ^n+i) = (c»+i/a»+ j) tanh (84), 

or, more conveniently, 

sinh (2i/'„+, - = (bn/a„) sinh i/^n, cosh (2i/'„+i - = A„/a„ . . . (85). 

As in Section VII (i), it is readily proved that 

^'Pot^o ~ dxpij^i = dxp^jS^ = ... = dxlf„lA„ = (86). 

If we construct the a.g.m. scale (ao=^ 1, h = k') and write 

... , tan f = sinh .Ao (87), 

it IS readily proved that 

Ao = A (i/r, A)/oos «/' and #o/Ao = #/A (i/-, k) (88). 

Thus, writing u = F(f, k), we have 

sn (m, k) = sin i/r = tanh i/r, , cn (u, k) = cos f= sech , 

dn («, k) = A (f, k) = A„/cosh x(f„ = cosh (2i/'i - i/'o)/cosh x[fj 
As n increases, the xff’s tend to a limit f„, while c„ tends to zero. 
Thus the result of integrating (86) gives 

f^n 

u = F(xf/,k) = (1/an) #»/cosh = (l/a„) tan"^ (sinh x/f„) . . .(90), 

If we notice that 

Aq dif/Q = A (i/f, A) ij/ 

= dn® (m, k) du/cn^ (u, A) = £* dn® (iu, k') du. 
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we have from 


(29) 

/Wo 

I A.^df^ = (\.li)Z{iu,U) + (E'IK')u. 


It is easily seen that equation (72) holds good when is written 
for Hence if we write 

'I'i) 


</ / / „ N _ o- tanh(2iA3-</^i) tanh(2iA:, 

2 *) cosh (2 i/'8 — >^ 2 ) *cosh(2i/'4-i^8) 


»*>■ 


and make use of (47) and (48), we obtain after some reductions 

Z(f,k) = tanh cosh (2\pi - - a„ sinh - 2' (fn, Cn) • • -(92). 

2 

The relations (27) and (28) enable us to calculate the second elliptic 
integral ^(^A, k). 


(ii) Calculation of the Jacobian Functions sn («, k\ cn (w, Jc\ dn {u, A), 
etc, in terms of the argument u. 

Compute the a.g.m. scale ao = h = determine from (90), 

which may be written 

sinh xf/n = tan (unu) (93), 

n being such that Cn is less than the small quantity determining the 
order of accuracy of the calculations. 

We then calculate the amplitudes ... \l ^29 from the recur- 

rence formula 

tanh (fn-i - = {On/cOn) tanh i/'n- 


Equation (89) then gives 

sn h) = tanh , cn {% k) = sech ^ 

dn {Uy k) = cosh (2xl/i - xl /^)/ cosh xj/^ j 

Equations (91) and (92) then enable us to calculate Z(Uy k). To 

compute <^{uy A), we note that (92) may also be written 

Z{Uy k) = sn {Uy k) dn (^4, A:)/cn (Uy k) - sinh 

~ i; (2^+^ - 1 ) (bn+i sinh - bn sinh xj/n). 

0 

Integrating this expression with respect to u and noticing that 

C^n 

sinh xj/ndi/JK = (l/^«) log cosh (2i/f„+a - i/^n+i), 

Jo 

we obtain the following series for @ (uy k), 

L® (0, i-COsh «/- J ^ ^ Lcosh (2xl/„+, - 

(95J. 


f 
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SECTION IX 


ON THE NUMERICAL COMPUTATION OF THE 
THIRD ELLIPTIC INTEGRAL 

As defined by Legendre, the third elliptic integral may be written 

n» = /o (1 Tn si 

The parameter n is taken as real, and may take any value, positive 
or negative. 

Writing u = jP(<^, h\n-- sn^ (a, k\ Jacobi considers the integral 

. f '^k^macnadnasn^udu . ,, r®(a-z4)“l 

n(u,a)= ^ 2 =t^^(a) + |log pr) ( ...(97). 

Legendre^s integral is connected with Jacobi's by the relation 

Tig (n, kf <l») = u -h sn a/(cn a dn a) . II (u, a) (98), 

and in this book is differentiated from it by the suffix 3. 

In order to include all real values of n in (96), a is not restricted to 
be real, and the treatment is usually divided into four cases : 

Case 1. negative between 0 and - 

Case II. negative between - 1 and - oo } Cases. 

Case III. n negative between and - 11 ^ 

Case IV. « positive between 0 and + oo } 

Case I. n negative, between 0 and - 1^. Hyperbolic case. 

We write ^ ~ F sn^ (a, ^) = — ^ sin^ 0, where 0 <a<K . . .(99). 

From (97) and (98) we have 
ww / 7 , N r.. sn (a, ^) ^ (a, k) 1 

- 1 1 (I'M)- 

^ cn (a, k) dn (a, k) ® L® (« + u, k)A ^ ' 

Compute the a.g.m. scale (a© = 1> bo=k'), and starting with <^o = make 
use of the recurrence formula (17), tan (<^w+i = (bnla<„) tan We 
then have sn (a, A) = sin <^oj on (a, k) = cos dn(a, k) = cos <^o/cos(</»i-<^o)> 

and Z’(a, A) =Cisin<^ + (?2sin<^2+ ••• +<?»8in<^„+ ..., while «na = <^»/2”. 
We thus have 
U^(n,k,<l>) 

tan cos ( ^o) ( ' j ’ j * I i”! 

1 + {^1 + ^2 Sin 1^2 + ... + Cn sin (f>n . . .jj 


Circular Cases. 


tan <^o cos (^1 sl^o) j 
* cos <f>o 
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The calculation of u by the recurrence formula (17), commencing with 
^0 = presents no difficulty. The most convenient method of computing 
the last term of (101) is to make use of Jacobies series 

® (iz?) = 1 - 2^^ cos 2ancc + 2q^ cos 4.ana! - 2(f cos + (102) 

in which we have written The nome q is most conveniently 

calculated from (64), log ^ = 2 log (ajan+i) - log (4ai/ci). 

If we write <l> - ^tt, the last term of (101) vanishes, and we obtain for 
the complete elliptic integral of the third kind the expression 

Ilg (», k, ^tt) = ir[l + tan ^0 sec <t>Q cos (</>i - <f>o) {cj sin + C 2 sin <^2 + 

... +Cnsin<3£»,,+ ...}] (103). 

We may also make use of the hyperbolic scale of decreasing amplitudes 
and the formulae of § 8. The final results are not, however, so con- 
venient for computation as (101) and (103). 


Case IL n negative, between ~ 1 and - oo . Hyperbolic case. 

We write 

w = — ^ sn^ (a + iK\ ^) = — 1/sn^ (a, A:) = — 1 /sin^ 0, where 0<a<K 

(104). 

Equation (98) becomes 

...( 105 ). 

We have 

^ (a + iK') = Z (a) 4* cn a dn a/sn a - ^ TrijK (106) 

and + iK') = iBH{z) (107), 

where log K' IK — ^irizIK, 

H{z) is Jacobi’s Theta-function given by the expansion* 

H{z) = 2^ (sin anZ-q^ sin sin 6anZ - q^^ sin lanZ + . . .) 

(108). 

Writing successively (a-^^) and (a^u) for 2 ; in (107), and a + iiT for a 
in (97), we obtain 

n (t^, a + iK') = uZ (a, A) + cn a dn a/sn a + ^ log [H(a - u)IH{a + v)\ 
and hence from (105), u<a, 

TL^in, k, <!>) = — sn a/(cn a dn a) [uZ (a, ^) + ^ log {H (a — u)IH {a + u)W 

(109). 

* Jacobi, Fundamenta Nova, §§ 62, 63 ; Ges, Werke, 1881, Bd. ii. pp. 229-231 ; 
Cayley, Elliptic Functions, 2nd Edition (1895), Chap. vi. 
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For numerical calculation, form the A.G.M. scale («o=l> and 

starting with 0o = make use of the recurrence formula (17), 
tan “ <#>n) = (hnldn) tan 

We then have ana = <f>j 2 ^y and as in Case I, we derive finally 

II 3 (Uf k, 

= ~^^tan 4 >q sec <#>0 cos {cisin + <^3 sin <^2 + ••• +Cn8in<^n+ •••} 

- 1- tan ^0 sec <#>o cos (<t>i - <^o) log [H(a-u)/H(a + w)] (HO). 

Starting with <^o = </> and making use of the same recurrence formula as 
above, the computation of u presents no difficulty. As before, the 
calculation of the last term of (110) is best achieved by computing the 
nome q by (64) and making use of (108). 

If we write = Jtt, the last term of (110) vanishes, and we obtain for 
the complete elliptic integral of the third kind the expression 

k, |7r) 

= -irtan^osec <^ocos (c^i - </io) {ci sin + C2 sin <^2 + +Cnsin<^n+ •••} 

......(111). 

When w>a, it will be noticed that the integral (105) has a singularity at 
«4=a, where the integrand becomes infinite. At the same time the logarithm 
of a negative quantity appears in (110), since H {a — u)-- H{u-a). 
It is necessary in this instance to define the path of integration, re- 
garding (105) as an algebraic integral obtained by writing z=^mu. Near 
the point z^sna this integral is approximately 

Ilg = sn^ a/(cn a dn a) dzl(sn^ a ~ z^), 
which integrated round a small semi-circle below the real axis gives rise 
to the imaginary term — ^rri m a/(cn a dn a). If we define the path of 
integration of the algebraic form of the integral as consisting of this 
small semi-circle and the remaining portions of the real axis, the above 
imaginary term cancels the term — sn a/(cn a dn a) . log (—1) appearing 
on the right-hand side of (110) when we write ll(a-u) = -l[(u-a). 
With this meaning attached to the path of integration, we may write 
the incomplete integral in the form, u> a, 

Ha (n, ^, <^) = - sn a/(cn a . dn a) [uZ (a, ^) -i- log {JS (u — a)/II (u + a)}] 

( 112 ). 

Numerical computation proceeds in the manner already described. 
Case III. n negative^ between - and - 1 . Gircula/r case. 

First Method. 

We write 

w = - P sn^ (ia + JT, k) = - (a, K) = - k%l - P sin® B) . ..(113), 

whence sin ^ = sn (a, k'\ 0 < a < K\ 
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Writing fa + if far a in (98), we find 

ITg (n, k, 4) = u -hi dn (a, ky[k '^ sn (a, k ') cn (a, h')\ . n (w, ia + K) 



while (97) gives 

n {u^ ia-h-K)=- uZ (ia + JT) + 1 log [® (ia - « + K)l® (ia -hu-h iT)] 



Making use of the addition formula (79) for Z (u), we find after some 
reductions 

iZ (ia + jS^ = iZ (ia, k) + ^8n (a, ^y[cn (a, k') dn (a, k')] ...(116), 
or, from (47), in terms of a real argument 

iZ(ia + JT) - Z(a, k') + ^ iraKKKJ ~ ^ sn (a, k^) cn (a, kyAn (a, k') 

( 11 '^)* 

Introducing Jacobies Theta-function ®i (z) defined by the relations 

®(z-h£:)=®,(z) (118), 

(5?) = 1 + 2^ cos + 2^ cos + 2^ cos + (119), 

we may write (114) in either of the forms 

^ . N dn (a, k') 

n,(», *, ^)=FiE^).cnCa,>l') 

^ / T . N dn (a, M) 

For numerical computation we may make use of two sets of formulae, 
according as we employ circular or hyperbolic functions. 

(i) Use of circular recurrence formulae. 

Form the a.g.m. scale («» = 1, Commencing with < p „ = < l >, make 

use of the recurrence formula (17), tan (<^„+i - ffin) = (bjon) tan 0„, and 
thus derive 

Commencing with i^o = ^, employ the recurrence formula (67), 
sin (2^,+, - ^f-,) = ( V«») sin 1^*, gi’'dng a„a = logtan (Jw + The 

application of (47) and (76) to both (120) and (121) gives us the formula 

^ f ^ JS^)] 

( 122 ) 

for the numerical calculation of the incomplete integral, and the formula 
n, (», k, Itt) = 2 (JT/P) cos (2i^i - V'o)/8m 2i/', . [a„ sin i/', + S (i/')] 
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The peculiar difficulty arising in the numerical computation of the 
circular cases of the incomplete elliptic integral of the third kind lies in 
the evaluation of the Theta-functions of complex argument. Writing 
= M + in Jacobi’s expansion (119), the series may be resolved into 
real and imaginary parts*. If we write 

tanS^- 2anU . sinh sin ^anU . sinh 4awa + . . . 

1 + 2g cos 2ant4 . cosh + 2g^ cos . cosh + . . • 

(124), 

we easily find log[©i(M-za)/®i(t4 + ia)] = 2e'^, from which it follows 
that (122) may be written 

Ilg (w, Ic, (fi) = (2/^'^) cos (2i/^i - »/<o)/sin 2i/^o b i^n sin + 2 («A)} - 

• (125). 

Powers of the nome q arising in (124) are most conveniently calculated 
from (64), 

to 

log ^ = Z (^)’‘“‘' log (d^K+i) - log (4ai/ca). 

The use of hyperbolic functions in (124) may be avoided by noting 
that if we write 

2cosh2a,ta = Q^+ 2cosh4an0t= + . 

2 sinh 2(3^»a = ’ 2 sinh 4ana = ‘ ® 

from which it follows that (124) may be written 

.T. _ ^(^^ ~ 2^ + - Q~^) sin 4^ + ~ Q~^) sin 6^ + . . . 

~ 1 +g(Q^+Q“^)cos24>+^(Q^-HQ“^)cos40+g^®(^+Q''*)co864»+... 

(124 bis). 


(ii) Use of hyperbolic recurrence formulae. 

Write n = — k^s>T)^ (^a + Kyk) = — ]^l{\ — h'^ tanh^ <^o0 (126)j 

whence sn (a, Ic) = tanh 0 < a < K\ 

Prom the hyperbolic recurrence formula (37) 

tanh (<^'n+i - <i>n) = (&w/»n) tanh <i>n 

we have 

ce,,a = ^„ 72 '' = sn (a, A') = tanh 
cn (a, ^') = sech , dn (a, K) = sech (</>i' - c^©'). 

* This procedure is due to Legendre, Fonctions Elliptiquesy Deuzi^me Supplement 
(1828), § IX. p. 146 et seq. 
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Thus from (44) and (120) we obtain 

TT (- I. J.N _ 1 

’ W - ^'«tanh <!>; cosh(^i' - <^o') 

r Msinh^o' ^ • v I ^ 1 •! ®i(M-*a)'] /,„„s 

To calculate u we make use of the recurrence formula (17) 

tan - <l>n) = (bjan) tan 0^, giving a^u = ^^2”- = <E>. 
Equation (124) then becomes 

,p __ 2g^ sin 2^ sinh 2^' + 2^ sin 4^ sinh W + 2/ sin 6<^ sinh + . . . 
1 + 2^ cos 2<E>cosh 2^' + 2g^cos 4^ cosh 44>' + 2/ cos 6^ cosh 6^' + . - 

(128), 

Equation (127) then gives 

nr h rfr\- ^ cosh 

s ) 9) ^/2 



and 

n, (w, h\^)=Y^ tanh^o'^co^W^) 

(129 *^5). 

Second Method, 

If we write P-K' -a, (113) takes the form 
71 = ~ ^ sn^ {K + iK' - ip) = - dn^ (yS, ^') = - 1 -f sin^ 0</3 

(130). 

We easily find 

= M+ MM) 

sn (/?, i:') cn (/8, ^') 

.[»z{i-*.x-i«-5.+i;iog (131) 

dn()8,^') 

~ k'^ an 

- [- & " i- «■ 

where -£ri(;?) is Jacobies Theta-function defined by the series 

111 (z) = 2gi (cos anZ + cos 3a«;2? + cos 5anZ + cos -t- . ..) 

(133). 
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Writing log [Hi (u + (u - (134), 

we find 

tan + ^mii5anU&mYi5an^+ ♦ - 

~ cosanUCOshanP^q^oo^SanUCOshSanfi+^oosbanUCOshoanP’^--^ 

(135). 

The numerical computation may be carried out by two methods. 

(i) Use of circular recurrence formulae. 

Starting with \l/o = 0' in (130), make use of the recurrence formula (67), 
sin (2il/n+i - if/n) = {hnjan) sin if/n- From (80) we have 

iZ (K + iK' - ip) = an (1 - sin lAn) ~ 2 (i/^), 

while 

sn(/3, ^') = sini/'o, cn (/?, ^') = cos > dn (^, ^') = cos ( 21/^1 - 

Thus (131) becomes 

Hj (», k,<t>) = u + [^' - « {«« sin ^„ + 2 (>/')}] ...(136). 

As before, u is most conveniently calculated from the recurrence 
formula tan (<l>n+i - cl^n) = (bnM tan <Jbn, giving anU = <#>n/2’‘ = 

The use of h 3 rperbolic functions in (135) may be avoided by writing 
log Q' = an ^ = log tan (irr+ |i/^n), so that 

2coshan^= ^' 4 - 2 cosh 3an^ = Q'® + Q'-'l . 

2 sinh anP =Q' - * 2 sinh ZanP = Q'® - * 

Thus (135) may be written % 

^ {q - q-^) sin (Q '^ - Q'-'^) sin 3 ^ + ^ Q -^) sin 5^+ . . . 

(Q' + Q'-^) cos 4> + (Q'^ + cos 3<1> + ^ (Q'* + Q'-^cos 5^+ . . . 

(138). 

Writing u = Kiji (136), and we have, since ^7r = anK, 

n, (n, A, Jtt) = iT [1 + K (1 - sin fn) - 2 (^)}] .. .(139). 

(ii) Use of hyperbolic recurrence formulae. 

Write sin B* = sn (/?, K) = tanh Q'lid (130) becomes 

7 ^ = - ^ sn^ (iST + iK' ~ ;/3) = - 1 + tanh^ <^o' • • .(140). 

Make use of the hyperbolic recurrence formula (37) 
tanh (^ n+l “ ) ~ (pnlar^ tanh <f>n } 

and we have 

an^ = <>' = <l>n/2\ sn (P k') = tanhe^o', 

cn (P M) = sech <^o', dn ()8, ^') = sech (<^' - <^o')* 

Also - iZ (iPi h) - Cl sinh sinh <^ 2 ' + . . . + Cn sinh <f>n-¥ .-.y 
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SO that (132) gives 

Hs (w, h, <l}) 

1 cosh 4*0 


cosh {4*1 -4^0) tanh 4>o 


ui^ sinh 4>o cosh {4>i— 4>o) + S sinh<^„' + . . .(141). 

As before, we calculate u from the recurrence formula 
tan {4>n^i - 4>n) = {bnM tan 4>n, 
giving anU = 4>n/^”' = Equation (135) gives for 

tan ^ ^ + sinh 34»' sin 3^ + sinh 5^^ sin 5<^ + . . . 

cosh 4>'cos 4> + (/^cosh 3^' cos 34» + o®cosh 54>' cos 54» + . . . 

(142). 

Writing u = Km (141), we have ^' = Jtt, and since = 

IIs (7^, ky ^Tt) 

(H3)- 

Case IV. n positive, between 0 and x . Circular case. 

First Method, 

Write F sn^ (^a, K)-W sn® (a, A')/cn’^ (a, ^') = W tan^ B,,, (144), 
giving sin ^ = sn (a, 1^\ 0 < a < K\ 

Thus 

XX f 1 JL\ sn (a, ^') cn (a, ^') f . V 1 ., ®(u-ia)l 

n. i, ♦) . » + - Vsi;) r 2(») + !• >»« 

(145). 

@(a?) is Jacobi’s Theta-function defined by the series 

® (z) = 1 - 2g cos 2anZ + 2^ cos AanZ - 2^® cos QanZ + . . . (146). 

Writing ^ log [0 {u - ia)/© {u + ia)] (147), 

we find 


ip _ 2g sin 2an u sinh 2an a — 2q* sin 4^^ ^ sinh 4a ,t a + ... 
1 - 2g cos 2anU cosh 2ana + 2^^ cos 4«„ w cosh 4a„a - . 


.(148). 


The numerical computation may be carried out by two methods. 

(i) Use of circular recurrence formulae. 

Make use of the recurrence formula sin(2i/^n+i~i/^n) = (5n/^*ji)sini/^n, 
commencing with 4^0 = ^ in (144). We then have 

sn (a, k') = sin ij/Q, cn (a, k') = cos 4/0, dn (a, k') = cos (2f 1 - 4/0), 
and by (78) 

iZ («a, k) = an. sin + 2 {4^) — sn (a, k') dn (a, ^')/cn (a, k') . . . (149). 

2 
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Hence (145) may be written 

( 150 )- 

As before, u is calculated from the recurrence formula 
tan - <t>„) = (6„/a„) tan </>„, 

giving (»).«< = <^»/2’‘=<t, while o„a = logtan(^ir + We have 
2co8h2a„o = ^ + Q-») 2cosh4a„o=Q* + Q-M , . 

2 sinh 2a„a= j ’ 2sinh4a„a= j ’ ^ 

where ^ = tan tt + 

so that (148) becomes 

tan u. - g ~ 2^ -q*(Q^- Q-*) sin 4^ + / (^ - ^-‘)sin6^- ■ 

1 ~3'(^+ Q"^)cos 2^+^(Q^+ Q-*)cos 44> - g® (Q* + Q"*) cosed*- . . . 

(152). 

Writing M = A" in (150), 4* = 0, and 

(153). 


(ii) Use of hyperbolic recurrence formulae. 

Write sin ^ = sn (a, ^') = tanh (1^^)> 

so that (^a, k)-1^ tan* 0 = A* sinh* 

Starting with <t>o given by (154), make use of the recurrence formula 
tanh (<^',H-i - ^b') = {bnjan) tanh <^*' (155). 

We then have 


sn (a, k') = tanh 

cn (a, k') = sech <^o', dn(a, k') = sech (<^/ - 


00 

and - iZ (ia, k) = ^Cn sinh 

Equation (145) may then be written 

Hs («, <A) = M + '^"co8h^^o^‘ ~ ~"~~ 

As before, u is calculated from the recurrence formula 
tan (<?!»«+i - 4>n) = (bjan) tan 

giving a^u ~ Equation (148) then gives for ^ 

2q sinh 24>' sin 24> - 2(f sinh 4^' sin 44» + 2j^ sinh 6^' sin 6<l>- ... 
^ 1—2^ cosh 2#' cos 2^+2^ cosh 4^' cos 4^— 2^®cosh 64>' cos 6<l> + . . . 


( 157 ). 
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Writing u = Km (156), ^ = 0, and 

n.(.. i,) = ir[i .(158). 


Second Method. 

If we write P = K'-am (145), we have 

« = -A“sn"(*Z’-*y8) = cn^(i8,/t')/sn'‘(/8,A') = cot“^ ...(159), 
giving m{S,k') = smB, 0<^<K'. 

We easily find 

Tl3(n,k,4>) 



where lT(z) is Jacobi’s Theta-function defined by 
ir(z) = 2q^ (sin a^z-f sin M^z + f sin har,z - (f‘ sin la^z + . . . ) 

If we write 


(161). 


we find 


= \i log \H{u + i^Wir u^i/3)] (162), 


cot = co^cCnUsinhanP - c os sinh +g^^C 08 5anU sinh . 

sin a^u cosh sin 3anU cosh Sanl^+ / sin 5an«^cosh5a„/8- . 

.(163). 


The numerical computation can be carried out by two methods. 


(i) Use qf circular recurrence formulae. 

Commencingwithi/'o=^in(159),make useof therecurrenceformula(67) 
sin (2«/^,,+i - = (5„/6in) sin 

From (78) we have 

- iZ {ifi, k) = tan \[/q cos (2fi-^o)- »n sin - S (i^) . . . ( 1 64), 

2 

while sn (^, k') = sin ifr^ , cn ((3, k') = cos ij/o , dn (/S, k') = cos (2»/^i - iI/q), 
Equation (160) then gives 

n8(w, k, <l>) 
siu 

^ 2cos(2f^-xl^,) ~ - S (tA)} + 

( 165 ), 

As before, u is calculated from the recurrence formula 
tan (^„,.i - <#>„) = (b„/a„) tan <Pn, 
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giving = O, while a„/3 = log tan (Jir + Thus 

2 cosh 3anj8 = Q'*+ 

2sinh3a„(8 = Q'®- 
Q' = tanQir + i«^„). 


2cosha„y3=Q' + Q'-M 2cosh3a,j8= 

= Q' - ^'■*1 ’ 2 sinh 3a„(8 = ^166), 


2 sinh a„/8 = Q' — Q'' 
where 

Thus (163) becomes 

tan cos 3^ + g® (C'" - Q'-')cos54'- . . . 

(O' + sin O - g2 (O'® + sin 3^ + g“ (Q''+ 4'-“) sin 5*-... 

Writing u-K'm (163), »' = ^ir, and (165) becomes ^ ^ 

Ha (w, ^^) = 2 CTO (2^ 7' -V 5 ~ a„(l-8in./'„)-S(^)] 

(168). 

(ii) Use of hyperbolic recurrence formulae. 

Write sin O' — sn (/S, ^') = tanh so that (159) gives 

n — —k^ sn^ (iK' - ip) = cot^ O' ~1 /sinh^ <^o', 0<P<K' ... (169). 

Commencing with </>o' in (169), use the recurrence formula 

tanh (<^ n+i - <f>n) = {bnjar) tanh </»n' (1 70). 


sn {P, k) ~ tanh <^o'> 
dn {p, k) = sech (<^/ - 


We then have 

cn (p, k) = sech <l>o', 

and -iZ (iPy k) — ^Cn sinh 0^'. 

1 

Equation (160) then gives 

n, (n, k, <!>) = 2 c„ sinh <f>n' + ^'~^ ... (1 71). 

As before, u is calculated from the recurrence formula 
tan (^„+i - <l>n) = (bjan) tan 

giving anU = <?5>n/2” = Equation (163) then gives for the expression 
cot ^ ^ — sinh S<^' cos 3<^ + sinh 5^' cos 5<1> - . . . 

~ cosh<^' sin ^ — q^ cosh 3^' sin 3^ + q^ cosh 5^' sin 5<> — . . . 

(172). 

Writing = A'in (163), ^' = Jtt, and (171) becomes 
ns («, k,^^) = K tanh .^o'^J^(^.^,'- .j.o;) ^ ^ 
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SECTION X 


NOTE ON THE CALCULATION OF THE THIRD ELLIPTIC 
INTEGRAL IN TERMS OF THE COMPLEMENTARY A.G.M. 
SCALE 


When k is very nearly unity, it may in some oases be advantageous 
to make use of the more convergent complementary scale («»'= 1, K=k'). 
The problem is to express the various functions appearing in the formula 

sn(a, [ rrf 7\ 11 ®(a-u,k)~\ 

I (a, dn (a, k) r ^ ^ ®(a + «,>l)J 

( 174 ) 


Ilg (w, A;, <^) = w + 


in terms of recurrence formulae involving the complementary a.g.m. 
scale. The four cases discussed in Section IX, depending on the values 
of n, require separate treatment. The reader will have no difl&culty in 
adapting the formulae of Section IX to the complementary scale should 
the occasion arise. We may note, in passing, the formula 


Q (g ~ cn (a + u) @ {i (a - u), k '} 

®(a + e^) cn(a — u)'&{i(a + u)yk'} 

which will be found to be useful. 


( 175 ) 



SUMMARY OF FORMULAE 


NUMERICAL CALCULATION OF ELLIPTIC FUNCTIONS 
AND INTEGRALS 

[Sections containing nm formulae are denoted by an asterisk *.] 
First Elliptic Integral. u=F(<l>,k) = 

Second Elliptic Integral. A(<^, k) = j* J(1 - k^ sin’’ E= E{\ir,k). 

Third Elliptic Integral. n,{n,kA )- [ (i ^ ^ sin^) V(1 -/^sin-^) ' 
Jacobian Functions. 

sin <^ = sn {u, k)y cos <^ = cn {% k\ sin^ <^) = dn {u, k). 

Scale of Arithmetico-Geometrical Means, 

«(, 1 h^ — k' Co -k 

«1 = i («0 + K) Cl = J (a, - 6o) 

«2 = J(«i + ^i) bi-fuibi C8 = J(ai-6i) 

a^i l^n Ofi 

1. Landen's Scale of Increasing Amplitydes 

tan (<^n+l - <^n) = {bnlctn) tan <l>n,. 

n>*>oo LJ^n ^ -J 

(K- E)IK= i {co^ + 2ci^ + 4c,* + ... + 2»c„* + ...), 

Z{d>,k)^E{<i>,k)-{EIK)P{<l>,k) 

= Cl sin 4>i + c, sin + . . . + c„ sin + — 
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*2. Calculatimi qf sn {u, k)y cn {u^ k\ dn (w, k) in terms of argument u. 
Calculate <!>„, - 2'^anU, n being chosen so that is negligible. Obtain 

••• <#»2, ^0 from recurrence formula sin <^^) = — sin</>n. 

an 

Then sn 1c) = sin cn (Uy k) - cos <^o) dn (w, k) = cos </>o sec (<#>i - <^o)« 


*3. The Hyperbolic Scale of Increasing Amplitudes 

tanh (<t>n+i - ^n) = (bjan) tanh <f>n^ Start with sinh </>o = tan <!>• 
F{<t>yK)= U [lK.<A„/2^, 


E{<i>, k!) = 


sinh <j!>o ^ 

cosh (<^i - <^o) 

“ (^1 


^F{4>,k’) 

sinh + Ca sinh </>2 4- 


c„smh<^„4- ...). 


*4. Calculation of sn (Uy k')y cn (w, k% dn (w, k') in terms of argument u, 

(1/i) sn {iuy k) = sn (w, A')/cn {Uy k')y cn (iUy k) = 1/cn {Uy k')y 
dn k) = dn (w, ^')/cn (u, k'). 

Calculate <l>n = ^ being chosen so that Cn is negligible. Obtain 

... <^o from recurrence formula 

sinh - <l>n) = (^?n/»u) sinh 

Then sn (t4, >t') = tanh <l>oy cn {u, k') = sech ^o, dn k') = sech (0i — ^o)- 

5 . Calculation of K\ E\ and the nome q— 

* Hyperbolic recurrence formula tanh Xn+i = V(^n+i/^n+i) tanh 2xn. 


k') = Lt [ 1 /^, 1 . tanh xq- Jk' sin </>. 

n-*.oo 


Put<^ = j7r. |7r^7ir=ilog(l/(?) = Ltr(i)Mog^1 . 

L t'w, Jn-SNoo 

(Legendre.) 

log 2 = 2 (!)”-> log ^ - log ^ . 

(Gauss.) 

iT' - £' = 1 i (1 - ... - -...) 


■<{lo8f‘-log|-llog|- 

...| -a„. 

6. Landeris Scale of Decreasing Amplitudes 
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*JS{>lr,k'} = anSmil>„ + ^(co^+2ei^+ ... +2®c„“+ ...) 


X — log tan (|ir + + S c„), 

2 


S ('/'». c«) = 2Ci 


tan (2i/ra - i/^t) ^ tan (21/^3 - }pi) 
cos (2 i^ 3 - l/fj) * cos (2 i^ 4 — l/fj) 


^ COS (2i/r„+3 - i/^„+2) 


^7. Calculation of sn A'), cn i'), dn (uy k‘) in terms of argument u. 

Calculate \l/n from sin = tanh n being such that c« is negligible. 

c 

Obtain ...fitiitfo from tan (^„_i - 1^.„) =-^ta,n>l/n. Then 

sn (Uy Ic) ~ sin cn (w, Id) = cos dn A') = cos (2i/^i - i/^o)* 

*8. Hyperbolic Scale of Decreasing Amplitudes 

sinh (2i/^n+i " W = ~ sinh iA». 

Ofyi 

F{\lfy &) = tan~^ (sinh tan ij/ = sinh 

E{^,k)-^F (i/', A) = tanh i/'o cosh (2fi - i/^o) - «» sinh - S' (i/f® , c»). 
Ji. 2 


*9. Calculation o/’sn (««, A), cn (t^, A), dn A). 

Calculate ij/n from sinh if/n = tan ^ being such that Cn is negligible. 
Obtain . .. ^A2, fi, from tanh = {cnjan) tanh Then 

sn {u, k) = tanh i/^o, cn {uy A) = sech ^^y dn {Uy k) = cosh - jf^o)/cosh 

*10. Calculation of Third Elliptic Integral — 4 cases. 


(W, 4 ^) - M + [wiT (a, A) + I log Q ^ ^j] . 


Case I. -k^<n<0. Case 11. — co < A^ < - 1. Hyperbolic Cases. 
Casein. -l<^^<-A®. r Method 1. (i) Circular, (ii) Hyperbolic 

recurrence formulae. 

Case IV. 0<w<+ao. "Method 2. (i) Circular, (ii) Hyperbolic 

recurrence formulae. 


Typical Solution, Case III, Method 1 (i). 

w = -A^/dn=*(a,A0==-A^/(l--A'^sin*(9), sin ^ = sn (a, A') 0<a<K\ 

1 . Form A.G.M. scale {a^-ly bo = k'). 

2. Derive UnU - </>„/2^ = ^ from tan (<^«+i ~ tan , commencing 

Un 

with <l>o = <l>- 
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3. Derive = log tan (Jtt + from sin (2\l/n+i - ^n) = ~ sin 

Ct/fi 

commencing with ij/^ = 0, 

4. Compute q from log g = S (J)”-' log («„/«„+,) - log (4ai/ci). 

5. Compute ^ from 

tan ^ 2anU sinh 2an€L + 2^ sin sinh + 2q^ + ... 

X + 2q cos 2anU cosh 2anOL + 2^ cos w cosh 4a^a + 2q^ + . . . * 

(Legendre.) 

6. Then, 

{n, k, i>) = gi„ + 2 

n. (», k, l^r) = + 2 



APPENDIX 


In the preceding sections only those formulae suitable for numerical 
computation have been given. Many interesting and elegant properties 
of the arithmetico-geometrical scales in conjunction with recurrence 
formulae not specially adapted to numerical calculation are given below 
as problems. In some cases the main steps of the proof are briefly in- 
dicated, and unless otherwise stated the final results are new. 
Properties of the a.g.m. scales (a© =1, = ^‘0- 

1. Prove that + Cn+i , bn = a»+i - Cn+i . Hence derive the series 

aQ = an + Ci + C2 + Cs + ..., 

bo=^an-Gi + C 2 + C 3 + Gauss, G, W. Bd. iii. p. 376. 

2. Prove that 

\/ Cn+2 ~ ^ ~ bn)j fj <35, ^+2 “ ^ ^ bn\ >/ Cn^2 ~ ^ n/ ^ n+a 

and hence \lan= — sjc 2 - s/c^— 

sjan^ + Jc 2 - .... Gauss, G. W. Bd. m. p. 376. 

3. Prove that + 2c^n+i = 2(an + c^+i, and hence derive the 

expansion 

{K - E)IK = (^0 + ai) Cl 4- 4 (^2 + ^3) C3 

+ 16 (<354 + af) C5 + ... + 2^” (a^n + Ggn+i) G2n+1 + • • . . 

4. Prove that 2bn-\Cn = — 2c\, and hence derive the series 

Co^ = 4 (b^Ci + 4^iC2 + 16^2^3 + . . . + ^'^^bnCn+i +.••)> 

an =60^ + 2 (60C1 + &1C2 + ••• + bnCn-hl + • • •)• 

5. From the results of Ex. 4, prove that 

EIK=bo + 2 ( 60 C 1 + 2 Z>iC 2 + 4^2C3 + ... + 2'^~^bn^iCn + ...), 

EjK =an^^2 (biC2 + + . . . + (2^^^ - 1 ) bn-iCn +...), 

and hence the identity 

l=an^ + 22 (2”-^- 1 i:2"Cnl 

2 0 

6. Prom the recurrence formula of Ex. 4 prove that 

2 (boCi + biC2 + ...) = Jco" - 1 (Ci^ + Ca^ + + ...), 

and hence = bo^ + 1 (ci^ + Ca^ + 6*3^ + . . . ), 

= Go" - - f \ci 4- C2^ + +...). 

Gauss, G. W. Bd. iii. p. 376. 

7. Derive the recurrence formula 

log (an-ijaf) = log (bnlhn-i) - 2 log {bn+ijbn) 



36 


APPENDIX 


and hence prove that 

\TrKIK^ log (4ao/co) - log (ao/^i) - i log (oi/a^) - . . . 

-(1/2^) log KK+i)-..., 
jK— log (4fei/Co) + ^ log (^2/^1) + I log (^3/^2) + • • • 

+ (3/2^) log (bnjb^+ .... 
Gauss, (?. W. Bd. iii. p. 377. 

[Note that ^ 7 rK'IK=(ll 2 ”') log (ajcn) as n increases indefinitely.] 

8. Derive the recurrence formula 

(1/2^-^ log = (3/2^-^0 log (ajbn) 

+ (l/2^-'') log (an/br,) - (1/2”) log (an^ilbn-i). 

Hence prove that 

= log (16ao VO + 32 (l/2”)log (^6,,). 

Jacobi, F. F. § 52, p. 201. 

9. Show that the a.g.m. scale (ao = 1, bo = k') may be calculated from 
the trigonometrical recurrence formula 

sin ^«,+i = tan^ ^ , 

commencing with sin ^0 = sin = (l - A')/(l + etc. 

Hence derive the following formulae : 

(i) Ci=^A;tan^^o) C2=j^tanA^otan|-^i, (?3=|A;tan^^otan|^itan|^2> 
c^ = sin , C2^ = sin 0 ^ sin 0 ^, ci = sin 0 ^ sin 6 ^ sin ^3, 

from which, in general, expressions for Cn and c^ may be formed. 

[Note that Cnjon = sin 0 ^^ bjan = cos 0 ^.'] 

(ii) log (2Z/7r) = 2 log (1 + sin = 2 2 log sec 

(hi) log {irl{ 2 Kk')} - J 2 log sec 
1 

(iv) \ttK’\K=- - i log g = -I log (4/sin ^1) - 2 (if log (1 + sin ^,,+1). 

[To derive (iv), make use of formula (61).] 

10. Show that the a.g.m. scale (ao =1, bo = k'^ may be calculated from 
the hyperbolic scale of decreasing amplitudes 

tanh 2xri+i = tanh^ Xn. or x«+i = \ log cosh 2xw, 
commencing with tanh 2xo = ^, or Xo = i log {(1 + k)l(l - yl;)}. 

Hence derive the following formulae : 

(1) tanh Xo , c^-^k tanh Xo tanh Xi , 03=^^ tanh xo tanh Xi tanh X2 , 

c-^ = tanh 2xi, ci = tanh 2xi tanh 2x2, 
tanh 2xi tanh 2x2 tanh 2xs, 
with similar expressions for Cn and V. 

[Note that b^ja^ - sech 2x71? = tanh 2x« •] 
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(ii) log(2Z/;r) = 2log(l + tanh2x«). 

(iii) logW(2ir-?;')} = 2ix». 

1 

(iv) \ttK'IK= log = log (2 coth Xo) - 2) ©"logCl + tanhx«+i). 

1 

[To derive (iv) make use of formula (61).] 

11. Show that the a.g.m. scale {a^ = 1, = Jc) may be calculated from 

the hyperbolic scale of increasing amplitudes 

taiih2i/^n = tanh^i/^^+i, or = |- log cosh 
commencing with tanh 2i/^o = or = i log {(1 + ^')/(l ~ ^0}- 
Hence derive the following formulae : 

(i) = i b,,la^ = tanh 2</^^ , cjan = sech 2\ffn . 

(ii) log (2jff7^) = 2 log (1 + sech 2i/^,,)* 

1 

(iii) log {7r/(2Z'A')} = i; log coth 

0 

(iv) ^kK'IK = log (f^- \l/n/2^~'^, as n tends to infinite values. 

[To derive (iv) make use of formula (60).] 

12. Writing 4 = {1 + (1 - show that the recurrence formula 

of Example 1 1 may be written 

2if/n+i - 4.ll/n = log 2 + log 4. 

Hence, denoting /3 = (1 - + Jk'), derive the expansions 

(i) log = log (2/iS) + 2 (!)”'■' log 4+1 . 

(ii) log q = log ap ) + 2 a^)^ -f 13 apT + (^pr + 

(iii) q =lP-^2 apy + 15 + 150 (1^)^^ + . . .. 

(See also Example 24.) 

13. Derive the following series for tt from the formulae of Example 12, 

TT = log^~^ — 2af^ — where ^ = |- (2^ - l)/(2i + 1). 

[Write k = k' = 1/V2, then K^K and log q~^ = tt. Note that 

= 23*14086, log ^-' = 3*141602..., 2^??^ = *0000070, 
so that two terms of the above series give the value of tt correctly to 
the tenth decimal. The remarkable approximation represented by the 
first term of the series was first pointed out by Legendre, F, E. t. i.] 

14. Given the nome q, compute the moduli k, k'. 

[The recurrence formula of Example 11 may be written 

4iA, - 2i/.,+i = - log 2 + log (1 - 

= — log 2 — — ... 
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Note that when n is sufficiently large = i log We then 

have 

Aif/n = 2” log q~^ - log 2 - 

Choose n so that the last term is negligible compared to the first two. 
The repeated application of the above recurrence formula enables us 
to calculate successively i/'n-i ••• »A 2 and finally xj/o, in terms of which 
k' = tanh k = sech 2i/^o*] 

15. Derive the differential relations 

^ .7-" -u«) = I- 

for all values of w. Gauss, W. Bd. iii. p. 380. 

16. From the results of Example 15 obtain directly the formula 

91-^00 f-'n/ 

[In terms of the a.g.m. scale (a^, bo, Cq) and the complementary scale 
(^o> formulae of Example 15 give 

But ao = ao', bo = Co', Co=^b(,' and a„^=bo'‘ + Co\ from which it is easily 
proved that 

Hence we conclude that d (~ log 0 + log = 0. 

We also have the identity d = 0, from which it 

follows that Cn~,, Cn being a constant independent of 

(«n, o„), {an, bn, Cn\ lu terms of Gauss’ notation for the a.g.m., 

we note from Section III that 

dn M {do, Co') — 2 M {dm Cn) — 2^ dnM ^ (ij Cnfdn^, 

while an = M{do, b^ = M{dm bn). Thus when n becomes indefinitely 

Jf(l, c„/a«) = 4- 7r/^ = |7r/log (4/^') as h' becomes small. Thus, in the 
limit as n increases indefinitely, C 7 i = -|7r, so that 
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17. Making use of the result of Example 16, prove that each of the 
ratios in Example 15 is given by 

V == Lt [^7rl(anan).dlog(anlan)]. 

n-^aa 

18. Derive Legendre’s relation EjK + E'jK’ — 1 = ^Trj^KK') directly 

from the properties of the a.g.m. scales («o= L h^k') and 1, k). 

[Differentiating the identity log (a^bn) - log = log(an/^«) and 

making use of the results of Example 15, we have 
d log (anbn) - d log = d log (an/bn) = 2^ (crt^/co^) d log (ao/bo). 

Writing down similar equations for ^^~l, n-2, ... 3, 2, 1 in place of 
n and adding, we find 

d log (a^bn) - d log (a^ b^) = (2ci^ + 4^2^ + . . . + 2^c J) Co"® d log (go/^o) (i) 
and 

d log (an bn) - d log (gq b^) = (2ci'^ + + . . . + 2'' c«'^) c^-^d log (aojbQ) 

(ii)* 

Since Go = Go' = 1 and bo = Co, we have 

d log (aobo) ~ d log (G o^o) = d log (Co/bo) = Gq^Co"® d log (go/6o). 

[Ex. 15 with 7^-0.] 

When n becomes large, an-bn^ so that 

log(G„&n)-ci^log (dnbn) = 2^1og (g^K') = - (ijir) G^^G^Co"^ C?log (Go/Jo)- 

[Ex. 17.] 

Note also the identity Co"^ <?log (go/6o) + Co'“^ log (aojbo) = 0. [Ex. 16.] 
Subtracting (ii) from (i) and dividing out by 

Co“^ d log (aojbo) = - d log (aojbo) 

we find 

(4/7r)G,,G,,'-Go^ 

- - (2Ci" + 4 c2^ + . . . + + . . . + 2^ + . . . ) 

(iii)* 

But we have (K- E)jK=^(co +2ci ... + 2^c„^+ ...), Gn= \^jKy Go= 1 
with similar expressions for the complementary functions in terms of 
accented symbols. Thus (iii) becomes 

Trj(KK’) - 1 = - 2 (1 - EjK) + Co^ - 2 (1 - KjK) + C 
or, noting that Co^ + Co'^ - 1, we derive Legendre’s relation in the familiar 
form 

EjK^ EjK - 1 = \itI(KK). 

The equivalent of the series (iii) is given by Gauss (G. W, Bd. ni. 
p. 380), but it does not seem to have been hitherto pointed out that 
Legendre’s relation follows immediately from it.] 
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19. Verify the following identities : 

(l + 2a^ + 2a^ + 2af + ...)(l-2a? + 2a^-2a^+ ...) 

(l+2a^ + 2iv‘^+2w^+ ...y + (l-2a^ + 2a}^-2af^+ ...y 

= 2(1 + 2aP + 2^ + 2cc^^ + • • •)^ 

(1 + 2a; + 2a^^ + 2a)^ + .,.y-(l-2x + 2^^-2^+ ...^ 

= (2a;i + 2a;^ + 2a;^ +.../. 
Gauss, Nachlass, G. W. Bd. iii. p. 447. 

20. Prom the results of Example 19, prove that if we form the 
A.G.M. scale commencing with 

ao = A (1 + 2^ + 2^^ + . . .y, h^ = h(l-2x + 2x^- . 

CQ~h {2x^ + 2x^ + 2x^' + . . . y^ 
the nth. term is given by 

an = h(l + 2Xn + 2^/ + ...y, hn=^h{l- 2Xn + 2Xn^ - ...)^ 

Cn^h (2Xn^ + 2^, ^ + 2Xn^~ + • • 0^ 

where log Xn - 2^ log x and h~M{a^, h^). 

Gauss, Nachlass^ G. W. Bd. iii. p. 448. 

21. Prove that x as defined in Example 19 is identical with the 
nome q when (ao= 1, io = ^0* 

[When n becomes large, 

-\\ogq = \irK'IK= Lt ^,log(a„/c„) 

n-^-Qo ^ 

= ^ log [(1 - + . . .yiiwj (1 + + • • OT 

n,-*.oo ^ 

= -1^1 ^ilogiT,. 
n-»*oo ^ 

= - J log iV, 

so that x=^q.\ 

22. From the results of Examples 20 and 21, derive the following 
expressions : 

/p = 1 -2^ + 2(y^-2g^+ ... . , __ 2q^ + 2q^ + 2q^^ + . . . 

1 + 2q + 2q^ + 2g® + . . . ’ ~ 1 + 2q + 2q^ + 2q^ + ’ 

J7r(l + + 2g^ + ...)^ KK = |7r(l - 2g + 2g^ - ...)l 

[From Example 20 we obviously have c»o= 1, Co = ^, ^ = 

23. If we denote ^ = (1 - n/^')/( 1 + V^'), derive the following ex- 
pressions : 

= {q + + . ..)/(! + 2q^ + 2g^® + ...), 

.r = 27r (1 + 2g^ + 2g^« + ...)V(l + sJWy = Stt (g + g® + ^ 
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24. From the preceding examples derive the well-known expansion 

+ 2 + 15 (J/3)® + 150 + . . . etc. (See also Example 12.) 

25. If we write 

a„= 1 + 2^n + 2^/+ ..., )8n = l-2^« + 2^„-*- 

yn = + 2Xr^ + 2Xn^' + ... , 

where log = 2^ log q, prove that the recurrence formula 

, ^ , tani/^^+i = (a^+,/^^+i)tan2«/^^ 

leads to the result 

Lt (^nlr)=^(^7r/K)F(i>,kX 

n~^ao 

where tan tan and K is given by the ^-series of Example 22. 

26. In terms of Landen's scale of increasing amplitudes (Section IV) 
tan (<l>n+i - ^n) = (bjan) tan </)^, derive the following expansions : 

dn(w, A) = A(<^o,^) = «n + ^JiCos<^i + C2COS02+ ... + c^cos(/>^+ 

00 

^ sn umu = J sin cos ~ <#>^). 

27. When (Cn+i/an+O is sufficiently small, prove that of Example 
26 may be computed from the expansion 

^n+l — 2^,1, — sin 2<l>n 

+ i" (^n+i/^«+i) 40% — ^ (pn+it ^n+lY ^0n 

28. Prove that the recurrence formula of Example 26 is equivalent 

to = 1 (^\-i + and hence show that 

log [® (m)/® ( 0)] = 2 log (a„/A„). 

0 ^ 

[Jacobi, Fundamenta Nova (1881), § 52, p. 204.] 

29. Prove that the recurrence formula a% tan 4>n+i = ^n+i tan 0%, where 
A% = sJ{anQo^^ 0% + J%^sin^ 0%), may also be written in the forms 

sin 0%4.i _ tan 0%4.i __ s! a n+i 'Jbn + _ Jbn + A ^ 

sin0% a% + A,/ tan0% >/a% + A,i’ \/a% + A,j* 

30. Prove that as a result of the recurrence formula of Example 29, 
<^^n/^n = d4>n+il^n+i- Hcncc derive the following formulae : 

F{^yF) = Lt (0%/«n), where 0 = 0o, sin 0o = sn (w, k\ 

Z (< t >, k ) = 2ci cos ^0 sin + dca cos sin <^2 + . . . + 2’‘c„ cos^„_i sin , 
log [® (u, k)l® (0, A)] = 2 2” log [2a„/(a„ + A„)]. 
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[The recurrence formula which lies at the basis of the results obtained 
in Examples 29 and 30 was originally given by Gauss in 1818 in his 
paper ‘‘ Determinatio Attractionis, etc.” G. W. Bd. iii. (1866), p. 353, 
and was further discussed by Jacobi, “Zur Numerischen Berechnung 
der EUiptischen Functionen,'’ G. Bd. i. (1881), p. 351. It is easily 
seen that these formulae are not adapted to rapid numerical calculation.] 

31. If sin = sn (m, k) prove that 

A]^ A2 . . . A,^ • • • I I 

tan anU = tan . 

ao«ia2 

[The result follows from Gauss’ recurrence formula of Example 29 ; 
see also Jacobi, F, N,, §36, G, W, Bd. i. p. 154.] 

32. Prove that Landen’s recurrence formula 

tan (ij/n - iA»+i) = (Cn+ilan^i) tan 
may be written tan 2 x ^+1 = sjicnja^ tan Xn if we write 
tan ij/n = sJianjCn) tan Xn- 

33. If we make use of Landen’s scale of decreasing amplitudes 

sin - xpn) = {bnja^ sin ij/n, 

commencing with = i’r, prove that log tan (Jtt + when 

n increases indefinitely. Hence show that the nome q may be computed 

from the formula q = cot^ (Jtt + ^ij/n) = (1 -- sin tAw)/(l + sin ij/n). Show, 

furthermore, that Legendre’s relation leads to the identity 

an (1 - sin ij/n) = (W-^) {1 - tanh {irK'l2K)} = S c^), 

2 

where 2 (i/'n, Cn) is the series of equation (75). 

2 


34. Prove that if the a.g.m. scale (^o, is computed, the circum- 
ference I of an ellipse of semi-axes is given by 

I = (2Tr I an) (ai^ — 2Gi - -...)• 


35. If we denote by and rg the greatest and least distances re- 
spectively of a point on one circle to the circumference of a coaxial one, 
show that in terms of the a.g.m. scale (»o = ^0 = ^ 2 ) the coefficient of 

mutual induction is given by 

Jf = {27r2/a,i} [ci + 2Gi + 4Cs® + . . . + 2^~^g^ +•••]> 
or, when the circles are very close together, in terms of the comple- 
mentary scale (ao' = n + ^ 2 , - rg), by the formula 


^ = 2TTan 


(r/--2^-4cs'^ -...). 


■ilog^,-ilog 



[King, L. V., ProG. Boy. Sog. A, Vol. c. 1921, p. 63 et seq. Several 
numerical illustrations are given in this paper.] 
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